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February 7, 1884. 

THE PRESIDENT in tlie Chair. 

Dr. David Gill was admitted into the Society. 

The Presents received were laid, on the table, and thanks ordered 
for them. 

The following Papers were read :— 

I. “ On the Motion of Fluid, part of which is moving Rota- 
tionally and part Lrotationally.” By M. J. M. Hill, M.A., 
Professor of Mathematics at the Mason Science College, 
Birmingham. Communicated by Professor G. G. Stokes. 
Sec. R„S. Received December 29, 1883. 

(Abstract.) 

Clebsch has shown that the components of the velocity of a fluid 
u, v , w, parallel to rectangular axes x, y , z, may always be expressed 

thus—■ 


u= f] *+\ d X v=^+\ d -±, w= C h+X d ±-, 

dx dx dy dy dz dz 

where X, Y are systems of surfaces whose intersections determine the 
vortex lines ; and the pressure satisfies an equation which is* equiva¬ 
lent to the following— 



where p is the pressure, p the density, and Y the potential of the 
forces acting on the liquid. 

It is shown in this paper that an equation in X only can be obtained 
in the following cases (that is to say, as in cases of irrotationa] motion, 
the determination of the motion depends on the solution of a single 
equation only) :— 

(1.) Plane motion, referred to rectangular co-ordinates x, y. 

d\ 

X being a function of x, y , t, let — be expressed as a function of 

dy 

X, x , t by substituting in it for y its value in terms of X, x , t found 


# British Association Report for 1881, p. 62. 
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from the equation which gives X as a function of x, y , t. Let the 

result be denoted by( —\. 

J KdyP 

Let differentials of X, x , t , regarded as independent variables, be SX, 
Sx, St respectively. 

Sx 


Thus 1 / dXV) is an integral taken with regard to x, X and t being 

w 

considered constant. 

S r Sx 


St I /is the differential coefficient of the above integral with 
J \dyf 

regard to t , X and x being considered constant. 

dX 
dt 

denotes that ^ is expressed as a function of X, x, t, and 
dy 

the result integrated with regard to x , X and t being considered 
constant. 

p and H are the symbols of arbitrary functions. 

Then the equation in X is shown to be— 

'dX V- 



d>\j8F(X,Q L, 
\ dx 3 dy^J j St 


dt 

dX 


( 


. =H 


X, 


f £F(X,£) j* Sx 

< sx~ 


(* Sx "] 
\dy) * 


and the current function is 


(* 

\Sx 

' dX\ 
dt 

J 

dX 

\ d v j 


St 


Lor a vortex of invariable form moving parallel to the axis of y 
dY 

with velocity — (where Y is a function of t only) the equation in X 
dt 

becomes — 


(cP d3 y>F(V) _ H 
Kdx% dy v St 


and the current function is 


< SK 


,t) .f I 1 

J /dXY > , 

\dy) J _ 


SF(X,t) ___^dY 
St " dt 

(2.) Plane motion, referred to polar co-ordinates r, 6. 
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With notation similar to that in the preceding case, the equation in 
X is— 


(»+l 


d 


\dr 2 


r dr r 2 dO 2 ) 


*Ky}+U 

st J 


/ dxY 
dt 
dX 

\de ) 


The current function is 

B¥(X 


dt 


— + 


=H 


f dX Y 
dt 
dX 
dB 


\ r®(V). 

< BX 


f rSr " 

j> 

KdeP 


For a vortex of invariable form rotating about the origin with 
angular velocity —, this equation becomes 

Civ 


( d 2 _j_l d ^ 1 d 2 '\ftiF(X,£) _ r 2 doj\ _ f fS¥(X,t) 

I 7i ~rr 1 Jri K J " 

i L 


_____ ___ 

\ dr 2 ~ V r dr * T r 2 dd 2 )\ Bt 2 dt) J BX 

and the current function is 


(* rBr q 


^F(X',^)_ r 2 dot 

Bt ' 2 W 


(3.) Motion symmetrical with regard to the axis of 0 in planes 
passing through it, referred to cylindric co-ordinates r, g. 

The equation in X is— 


r 2 \dr 2 r dr dz 2 



/dx\ 

n 

K aF( >>*) + |Vgr 

dt 

* ( 
/ 

/] « J 

1 dX j 


l 1 

\Tz\ 

J 


=H 


J B¥(X,t) 

\~sxr m 



The current function is 

m f 

Bt 


rBr 


fdxY 

dt 
dX 
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For a vortex of invariable form which moves parallel to the axis 
of z with velocity this equation becomes— 

l(d* _l d f * f£F(V)_ 1 

r 3 \e£r 3 r cfe 3 / ^ £X J ^dX^jj j ’ 

and the current function is 

aF(X^_r 3 dZ 
St 2 dt ' 

Suppose that in any of these cases any particular integral of the 
equation in X is taken. 

It is shown that the components of the velocity can be expressed in 
terms of X and differential coefficients of X, and that the current 
function is also known. 

In the case of a fluid, part of which is moving rotationally and 
part irrotationally, the boundary surface separating the rotationally 
moving fluid from that which is moving irrotationally contains the 
same vortex lines, and may be taken as the surface X=0. 

Now, if the integral taken of the equation in X do actually corre¬ 
spond to a case of fluid motion in which part of the fluid is moving 
rotationally and part irrotationally, the most obvious way to find the 
irrotational motion will be to find its current function from the con¬ 
ditions supplied by the fact that the components of the velocity are 
continuous at the surface X=0. If after taking any integral of the 
equation in X it be found theoretically impossible to determine the 
current function of an irrotational motion outside the surface X—0, 
which shall be continuous with tbe rotational motion inside it, then 
the integral in question does not correspond to such a case of fluid 
motion. 

In this method no assumption is made as to the distribution of the 
vortex lines (as in the method of Helmholtz) before commencing the 
determination of the irrotational motion. 

On the other hand, it is not a particular case of a method applicable 
to motion in space of three dimensions. 

But it can be shown that Clebsch’s forms for the components of 
the velocity do also lead to a method which is applicable to the deter¬ 
mination of the irrotational motion when the rotational motion in 
space of three dimensions is known. 

For the rotational motion being known, the components of the 
velocity are known for this part of the fluid. Let the components 
of the velocity be expressed in Clebsch's forms, so that x* \ Y r are 
known. 

Moreover, let the forms be so arranged that the surface separating 

vol. xxxvr. u 
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the rotationally moving fluid from that which is moving irrotationally 
is the surface \=0. 

Then at this surface the components of the velocity are —, —. 

^ 87 dx dy dz 

Now, obtain in any manner a velocity potential 0 for space outside 
X—0 continuous with % all over the surface X-=0. This is theoretically 
possible always. 

If the velocity potential so obtained make the velocity and pres¬ 
sure continuous all over the surface /V=0, then a possible case of 
motion will have been obtained. 

The conditions to be satisfied in order that the velocity may be 
continuous at the surface \=0 are that there 

dx dx dy dy dz dz 

In order that the pressure may be also continuous, it is further 

necessary that all over the surface \ = 0. 

dt dt 

The most obvious way of obtaining the velocity potential will be to 
apply Helmholtz’s method of finding the components of the velocity in 
terms of the supposed distribution of magnetic matter throughout the 
space occupied by the rotationally moving fluid. 

It must, however, be remembered, as is remarked by Mr. Hicks in 
his report to the British Association on “ Recent Progress in Hydro¬ 
dynamics,” Part I,* “ That the results refer to the cyclic motion of 
the fluid as determined by the supposed distribution of magnetic 
matter, and do not give the most general motion possible.” It 
appears also from Examples I and III of this paper that it is not possible 
to assume arbitrarily the distribution of vortex lines, even when it 
can be shown that the equations of motion are satisfied at all points 
where the fluid is moving rotationally, and then to proceed to cal¬ 
culate the irrotational motion by means of the supposed distribution 
of magnetic matter. For in these examples, values of the components of 
the velocity of a rotational motion, satisfying the equations of motion 
throughout a finite portion of the plane of x, y , are found. Thus the 
distribution of vortex lines, and, therefore, that of the supposed 
magnetic matter over a finite portion of the plane of a?, y is known. The 
surfaces that always contain the same vortex filaments are found. 
Inside one of these the supposed magnetic matter is distributed, the 
current function at external points is calculated by Helmholtz’s 
method, and it is shown that the velocity so found is not continuous 
with the velocity at the surface, which separates the rotationally 
moving liquid from that moving irrotationally. 

Another way (suggested by Clebsch’s forms) of obtaining the 
velocity potential will be as follows :— 


# Report for 1881, Part I, p. 64. 
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Calculate tlie quantity p — — JL 

H J ' 4 i\d& df dz*) 

Treating p as tlie density of a material distribution inside \=0, taking 
no account of the value of p outside the surface \=0, obtain the 
potential of this distribution. Let the potential inside \=0 be an( ^ 
outside let it be 0. 

X will, in general, differ from %; first, because x ma y contain many¬ 
valued terms, which may be denoted by P satisfying Lap]ace’s equa¬ 
tion; and, secondly, because %—P may be the potential of a distribu¬ 
tion of matter, part of which is outside \=0. 

Accordingly, it is necessary to examine whether it is possible to 
find many-valued terms P satisfying Laplace’s equation such that 
x' + P = X- 

Then 0 + P will be the velocity potential of the irrotational motion, 
provided that it give zero velocity at infinity. 

The few illustrations which follow are a first attempt to apply the 
above theory to particular cases. 

Example 1. x 9 y being rectangular co-ordinates in a plane, 

Y any function of the time t, Y its differential co¬ 
efficient with regard to the time. 

/, a , b , constants, 

u, v the components of the velocity parallel to the axes, 
it is shown that • v = Y —(/)?? satisfy the equa¬ 

tions of motion. 

The surfaces containing the same particles are the elliptic cylinders 

(/)(5+ Cj 9^) =COMtant ' 

For a finite portion of the plane of as, y outside the cylinder 
— 1, the current function of an irrotational motion 

a 3 b 2 

continuous with that inside the cylinder is— 

a v 4 . yz ___ 

(/) —^—l°g(v / tt 2 +e + + e) 


+ 2 / r 2 ( e + a2+& % .+ ffl2 + 62 \+\ 

a 3 -& 2 \a 2 -& 2 V 2 A 2 7 2 / 

where e, v are the elliptic co-ordinates satisfying the equations 

+ (y-Y) 8 -1 

a? -j-6 6 3 + e ’ a 2 -j-v b' 1 4- v 

— b 2 <e< oo , — a 3 <v< — b 2 , b < a. 


and 
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Example 2. Kirchhoff’s Elliptic Vortex Cylinder. 
x, f are rectangular co-ordinates referred to axes rotating with 
uniform angular velocity w about their origin. 

/, a, b are constants. 

u', v' the components of the velocity parallel to the rotating axes. 

_ )=const. in this case (- u , v are velocities parallel fixed 
\dx dijj 

rectangular axes in space). 

It is shown that %'==?/'( ; v'=—- x'l^L — w) satisfy the equa- 

\b 2 / \u 2 J 


tions of motion. 

The surfaces containing the same particles are the elliptic cylinders 


x 2 if 2 , 

+ = const. 

a 3 b 3 


If — (which is— —— — then the motion mav be supposed 

ah \ (ct + b) 2 J J Ll 

to be given by the above values of u\ v inside the cylinder 


7 /3 

?1_4-L_=rl 
a 3 b 2 


and therefore to be rotational; but outside this it may be supposed 
irrotational, and its velocity potential 


ISM 1 


u 3 + fr 3 + 2e _> 

W (a 2 e) (b 2 + e)' 


- g’afr sin' 


-l 


^ a 2 + e 


where e is the root of 


+ — = 1, which lies between — b 2 and oo , 

a 2 + e b 2 +e 


a being > b. 

Example 3. With the same values for if v ( as in the last example, 
except, however, the relation between / and w, it is shown that an 
irrotational motion continuous with the rotational motion inside the 

x ’ 2 n ' 2 

cylinder —— \ can exist between this and a confocal elliptic 
a 2 b 2 


smooth rigid cylinder surrounding it; provided that confocal elliptic 
cylinder be made to rotate with the same angular velocity w. 

The e of the confocal cylinder is— 


. W 2 

where — oo < < —— 

f ab 


,/a 3 + 6 2 _<«\ 
A aW f ) 
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The components of the velocity of this rotationally moving liquid 
are— 

and -*(¥-")• 

The components of the velocity of the irrotationally moving liquid 
are— 



Example 4. The vortex sheets are coaxial circular cylinders, and 
the motion is everywhere perpendicular to the radius vector from the 
axis, and a function of the distance. 

Thus radial velocity=0, velocity perpendicular radius vector= 

-np- ' ' 

This example is given merely to illustrate the expression of the com¬ 
ponents of the velocity in Clebsch’s forms. 

It is shown that 

V,= "(' + FF) e X F " (r)+ ^' (r) ) 

x =e[F(r)-fF(r)+ f IM } 

+ 1\ F(r)F'(r)-i(F(r))» + lF(r)F(r)- f^(F(0) 8 } 
l r Jr J 

If F(a) = 0, and the rotational motion be supposed confined to the 
interior of the cylinder r=a, then a suitable value of the velocity 
potential at external points is— 

0=0(-aF(a)) + t { -i(F( a ))2-J^(F'(«))2 } 

Example 5. r, z being cylindric co-ordinates. 

Z, an arbitrary function of the time ; Z, its differential coefficient 
with regard to the time. 

a, b, a positive constants. 

t the velocity from the axis of cylindric co-ordinates in the direction 
of r. 

w the velocity parallel the axis of cylindric co-ordinates. 

It is shown that t=2 ar(z — 7i)\ w~Z — — Z) 3 — 46(r 3 ~-« 3 ) 

satisfy the equations of motion. 
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One set of the surfaces which always contain the same vortex lines 
is given by ar 2 (z~Z) 2 A &(r 2 —a 2 ) 2 == const. 

If the constant be<6a 4 , the equation represents ring-shaped 
surfaces. 

The above expressions for the components of the velocity determine 
a possible case of rotational motion inside the hollow smooth rigid 
surface of annular form whose equation is ar 2 (z~Z) % -\-b(r 2 —« 2 ) 2 = 
const., and which moves parallel to the axis of z with velocity Z. 

The values of %, X, yfr of ClebscVs forms of expression for the com¬ 
ponents of the velocity are calculated, but the writer has not succeeded 
in finding an irrotational motion outside one of the above annular 
vortex sheets continuous with the rotational motion inside it. 

Could the solution be completed it would amount to a discussion 
of the motion of an annular vortex, all of whose parts are of finite 
dimensions. 


Appendix. 

The paper concludes with an Appendix, in which has been placed an 
account of Examples 1 and 2, the calculation of the potential of the 

elliptic cylinder —+ A=1, the density of which at the point x, y, is 

CL rJ 0 


i(L~ l 

'&Xb 2 d' 


1 \_®£_ 

*2 )/aP y*V. 

la 2 b 2 ) 


(Although this density is infinite along the axis of the cylinder, yet 
the total amount of matter in any elliptic cylinder surrounding the 

■>2 2 

axis whose equation is —-f—=const., however small the constant 
a 2 b' 2 

may be, vanishes. Hence it is not singular that the potential is 
finite.) 

The potential inside is— 


c a—b . cab . n a 

-- xy + —— sm" 1 -— 

2 a + fr J 2 / t 2 


v x 2 Au 2 


a 2 b 2 


The potential outside is- 


cab ( 


_ a 2 + 5 2 + 2e 

a v / (a2+e)( 6 *+"0" 


A , cab' . 


Vx 2 +y 2 


a 2 -f e b 3 -f- e 


where 



